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Abstract 

The variable-phase approach is apphed to scattering and bound states in an 
attractive Coulomb potential, statically screened by a two-dimensional (2D) 
electron gas. A 2D formulation of Levinson's theorem is used for bound-state 
counting and a hitherto undiscovered, simple relationship between the screen- 
ing length and the number of bound states is found. As the screening length 
is increased, sets of bound states with differing quantum numbers appear 
degenerately. 
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The quantum mechanics of low-dimensional systems has become a major research field 
with the advent of growth techniques for the realization of semiconductor quantum wells. 
Almost all of the computational techniques developed for three-dimensional (3D) problems 
have already been extended to lower dimensions. The variable-phase method for the cal- 
culation of the scattering phase shifts was introduced long ago by Morse and AUis and 
expanded by Calogero [@], Babikov |^ and others @] but has to our knowledge never been 
used for analysis of any realistic two-dimensional problem. In this Letter we use the variable- 
phase method to treat both scattering and bound states on the same footing for the Coulomb 
potential statically screened by 2D electron gas. This potential plays a central role in the 
physics of semiconductor heterostructures, one of the most rapidly growing fields in con- 
densed matter physics. Despite having been studied for approximately 30 years this 
potential exhibits some peculiar features which have never been noticed until now, e.g., as 
the screening length increases, sets of bound states with differing quantum numbers appear 
degenerately. These degeneracies are not, however, the same as the degeneracies in the un- 
screened spectrum. We also report a hitherto undiscovered, simple relationship connecting 
the number of bound states to the screening length. 

In what follows we consider a simple case in which a screened charge resides in the 
same plane as a screening 2D electron gas. This geometry is appropriate to the problem 
of screened excitons or impurities in a narrow quantum well. We assume that an 
attractive potential is created by a point charge e at the origin and use throughout this 
paper excitonic Rydberg units where length and energy are scaled, respectively, by the 
effective Bohr radius a* and Rydberg Ry* . In these units and geometry the Thomas- 
Fermi expression for the statically screened potential in the electron plane is 

where is the 2D screening wavenumber. Eq. (|l|) is the 2D analogue of the Yukawa 
potential. Taking the 2D Fourier transformation of Eq. (|l|) yields in real space 

. (2) 
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where Joiip) is the Bessel function, and p = (x^ + y'^Y^'^ is the in-plane distance from 
the origin. 

To describe the apphcation of the variable-phase method in 2D we consider a particle 
moving with energy E = k"^ in the potential V{p) which has radial symmetry. Since the 
potential is symmetric, we can separate variables in the expression for the wave function: 

oo 

^{p,^) = J2 cos{mip) , (3) 

m=0 

where m is the absolute value of the projection of the angular momentum onto the symmetry 
axis of the potential. At large distances from the scattering center, the radial function 
satisfies the free Bessel equation, whose general solution is 

Rm{p) = Am[Jm{kp) cos 6m - Nm{kp) sin Sm] 

I 2 \'/' 

( ^ j cos(A;p - (2m + l)7r/4 + b^) , (4) 

where is the scattering phase shift Jmikp) and Nmikp) are the Bessel and 

Neumann functions, respectively. Both total and transport cross sections in 2D can be 
expressed via the scattering phase shifts in a simple fashion 0. 

In the variable-phase method and 5m are considered not as constants but as 

functions of the distance p . The phase function 6m{p) is the phase shift produced by 
a potential cut-off at a distance p. Then the scattering phase shift can be obtained as a 
large distance limit of the phase function 6m{p) , which satisfies the following first-order, 
non-linear differential equation originating from the radial Schrodinger equation [Q: 

d TT 

^5m(p) = - -pV{p)[Jmikp) cos5m{p) - Nmikp) sinSmip)]^ (5) 
with the boundary condition 

SmiO) = . (6) 

Eq. (|^) ensures that the radial function does not diverge at p = . Physically one can view 
the phase function as measuring the retardation of the scattering wave function due to the 



potential. From this view it is clear that Eq. (j^) is the correct boundary condition since the 
potential can produce no retardation at the origin. The total scattering phase shift is given 
as an asymptote of 6m{p) '■ 

Sm = lim 5rnip) ■ (7) 

For numerical convenience, instead of the boundary condition Eq. (Q), the small-p expansion 
is used 

Up) - - I v{py'--Up' , p - . (8) 







It is also useful to rewrite Eq. (Q) in terms of special functions |^,|| as 

Vs{p) = -2|i-|g,[Ho(M-iVofep)]| , (9) 

where Ho(gsp) and No{qsp) are the Struve and Neumann functions, respectively. Series 
representations of these functions for small values of argument and asymptotic expansions 
for large values of argument allow the accurate calculation of K(p) for all values of g^p . 
Asymptotic expressions for the potential are 



Vsip) - -T^ji + E (-ir + OUp)-''''] } , g.P » 1 , (10) 



and 

Vsip) - -2[p-' + g.ln(Cg,p)]+Ofep) , g.p < 1 , (11) 



where C = |exp(7) , and 7 is Euler's constant. The asymptotic expressions Eqs. (|TO|,|TT]) 
show that the potential Vs{p) satisfies the conditions: 

/oo 
V{p') dp' ^ 0, p 00 , (12) 
p 

and 

pV(p) ^ , p ^ , (13) 
which are sufficient to allow application of the variable-phase method for this potential. 



Figure |l] shows the k dependence of the phase shifts Sm obtained by the numerical 
solution of Eq. with the initial condition Eq. (^) for the screened Coulomb potential 
Vs{p) . The distinctive feature of this plot to which we draw attention is that in the low- 
energy limit, k , the scattering phase shift is an integer number of vr : 

lim Sjn = i^T^ ■ (14) 



This behavior may be understood recalling Levinson's theorem ||T0| which connects the zero- 
energy scattering phase shift with the number of the bound states for non-relativistic parti- 
cles in 3D. Recently Levinson's theorem has been discussed for Dirac particles, multichannel 
scattering, multi-particle single-channel scattering, one-dimensional scattering systems, sys- 
tems with non-uniform effective mass, and even for time-periodic potentials [0. However 
its applicability to the 2D scattering problem has not been considered yet. We expect that 
this fundamental theorem holds also in 2D in the form of Eq. ([1^ , where v is the number 
of bound states for a given m. A simple proof for m 7^ follows from Eq. in the same 
fashion as Calogero's proof of Levinson's theorem in the 3D case 0. For m = the proof 
is more complicated due to the logarithmic divergence of the Neumann function, N^lkp) . 
A rigorous proof of Levinson's theorem in 2D, based on analytic properties of the scattering 
matrix |13|, is needed. We have, however, verified the validity of Eq. ([T^ for all m values 



in the analytically tractable case of a circular finite potential well 0]. 

Despite its appeal, Levinson's theorem has not been widely used to enumerate bound 
states since there exists an ambiguity in the usual definition of 6m, being defined only up 
to mod(7r) However using the variable phase approach avoids this problem since the 

phase function is uniquely defined by Eqs. for all p. 

In Fig. the number of bound states, obtained as the low-energy limit of the scattering 
phase shift in units of vr , is plotted as a function of the screening length = l/g^ for the 
potential K . As the screening length increases, the potential supports more bound states 
and these new bound states appear at critical values of the screening length indicated by 
the steps. One can see from the location of these steps that the u-th bound state for a given 



m appears at the critical screening length, given by a simple formula 

(2m + z/ - l)(2m + u) 



' c 



z/ = l,2, .... (15) 



This intriguingly simple relation has, to our knowledge, never been reported despite numer- 
ous calculations of the binding energy in the screened Coulomb potential, since conventional 
numerical methods for the binding energy calculation fail for extremely shallow energy levels. 
Note that for m = the first bound state appears immediately at (r^)^ = 0, corresponding 
to the fact that there is always at least one bound state in any symmetric 2D attractive 
potential. Using Eq. (p!5D we can simply evaluate how many bound states the 2D statically- 
screened Coulomb potential will support for any value of m. 

The Bargmann bound condition re-stated for two dimensions is mu < . This 



gives a gross over-estimation of the number of the bound states; e.g., for = 4, m = 1 
this implies only that there are less than four bound states whereas there is in fact only one. 
In this sense the Bargmann condition is of limited utility. 

For many applications (e.g., in a partition function calculation) the value of the binding 
energy, not just the number of bound states, is important. The variable-phase method 
provides an elegant and efficient solution of the eigenvalue problem as well. To approach 
this problem we recall that for the states with negative energy the wavenumber k is 
imaginary, k = in , and we introduce the function /^^(p, k) vanishing in the origin and 
satisfying a non-hnear equation H 



lm{l^p) COS Prn{p, «) + -Km{np) Sill Prn{p, 1^) 



TT 



(16) 



where Im{i^p) and Km{np) are the modified Bessel functions of the first and second 
kind, respectively. The functions Im{i^p) and Kmi^p) represent two linearly independent 
solutions of the free radial-wave Schrodinger equation for the negative value of energy, E = 
—K^ , and cot pm characterizes the weights of the diverging (/„) and converging {Km) 
solutions as p ^ oo . For the bound state, the diverging solution vanishes, implying the 
asymptotic condition 

6 



Hmip-^00,K^) = {U- 1/2)71, U = l,2, ... . (17) 

Here i' numerates the bound states for a given m and (z/ — 1) is the number of non-zero 
nodes of the radial wave function. For numerical solution of Eq. (|I6|) instead of the boundary 
condition /im(0, n) = an approximate initial condition (analogous to the condition Eq. (^ 
for the phase function 6m{p) ) is used. 

The bound-state energies Em,u versus the screening length for the potential Vg 
are plotted in Fig. ^ For clarity all the curves are normalized by the values of the energy 
for the unscreened potential [|l^, Em,v{qs = 0) = — (z/ + m — 1/2)^^ = — (n + l/2)~^ . The 
principal quantum number n is given by n = v + m — 1 . We note that the unscreened 
eigenstates with the same value oi m + v are degenerate. A consequence of the relationship 
in Eq. ([TSD is that when the bound states first appear (at = (r^)^ ), eigenstates with the 
same value of 2m -|- v are degenerate. With increasing this degeneracy is lifted and the 
degeneracy for the states with the same values of n appears as -h- 00 . 

In conclusion, we have used the variable-phase method to study scattering and bound 
states in two dimensions. The 2D analogy of Levinson's theorem is formulated and verified 
empirically. A hitherto undiscovered, simple relationship between the screening length and 
the number of bound states in the statically-screened Coulomb potential has been found. 
Challenging problems such as a general proof of Levinson's theorem in 2D and the analytic 
derivation of Eq. (pTsD are clearly the subject of further research. 

This work was supported by the Engineering and Physical Science Research Council and 
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FIG. 1. Scattering phase shifts versus in-plane wave vector k (in units of inverse Bohr radius 
1/a* ) for a 2D particle in an attractive Coulomb potential screened by a 2D electron gas, screening 
length rg = I/Qs = 5a* . Numbers show m values. 
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FIG. 2. Number of bound states, calculated as the low-energy limit of the phase shift, plot- 
ted versus the screening length (in the units of Bohr radius a* ). The solid line corre- 
sponds to m = ; dashed line, m = 1 ; dashed-dotted line, m = 2 . New steps appear at 
(r,)^ = {2m + u-l){2m + u)/2 , v = 1, 2, .... 
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FIG. 3. The bound-state energies Em,u plotted versus the screening length for dif- 
ferent m values. The curves are normalized by the energies for the unscreened potential, 
Ejn^,y{qs = 0) = —{m + u — l/2)~^ . Numbers show m and u values. 
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